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At finite temperature and chemical potential, the leading-order (hard-thermal-loop) contribu- 
tions to the gauge-boson propagator lead to momentum-dependent thermal masses for propagat- 
ing quasiparticles as well as dynamical screening and Landau damping effects. We compare the 
hard-thermal-loop propagator with the complete large- A*'/ gluon propagator, for which the usually 
subleading contributions, such as a finite width of quasiparticles, can be studied at nonperturba- 
tively large effective coupling. We also study quantitatively the effect of Priedel oscillations in 
low-temperature electrostatic screening. 
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^ ■ I. INTRODUCTION 

T— I ■ 

J> ' At finite temperature and chemical potential, propagators of elementary fields are modified importantly at soft 
, momentum scales. At leading order, scalar fields acquire simple tliermal mass terms, but fermions and gauge fields 
develop additional quasiparticle branches with compHcated momentum dependent thermal masses for propagating 
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£^ ] modes. Landau damping cuts, and, in the case of gauge fields, poles for imaginary wave vectors corresponding to 
dynamical screening Q, Q, |^ |j, |^ |^ . 
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While the physical singularities of propagators in a gauge theory are gauge-fixing independent 0, naive perturbation 
. theory leads to gauge-fixing dependent results beyond leading order and requires (at least) resummation of the full 
nonlinear and nonlocal hard-thermal- loop (HTL) effective action (see e.g. Ref. 9j for a recent review). It has also 
been found that thermodynamic quantities like the entropy of QCD can be understood well down to temperatures 
T > 3Tc with couplings as large as 5 ~ 2 by using HTL quasiparticles in <i>-derivable two-loop approximations 
(see also One would therefore want to understand the underlying quasi-particle picture also at larger couplings. 

^ ' The success of those resummation techniques has been established so far only by comparison to lattice data, therefore 
'Tt ' an independent test is desirable. One possibility to go beyond HTL is the large fiavor number (Nf) limit 
J> . in which thermodynamic quantities like pressure or entropy can be calculated to all orders in the effective coupling, 
kjp ' at next-to-leading order of the 1/iV/ expansion (which is of order N^)- In the large- A^/ limit, the thermodynamic 
■ potential exhibits a non-trivial (non-monotonic) behavior when going from weak to strong coupling. As will be shown 
d ' in a forthcoming paper the latter is due to a corresponding nontrivial behavior of the next-to-leading order 
asymptotic mass of the quarks. However, the contributions of the gauge bosons, which are needed only to leading 
order in the large- iV/ expansion, are of similar magnitude, and it is worth to study this propagator and its properties 
in detail, in particular how it compares with the corresponding HTL propagator. 

The large- TV/ gluon propagator is essentially abelian and one-loop, and some comparisons with the HTL propagator 
have been worked out previously in Ref. [16], however without discussing renormalization scale dependences that affect 
any approximate result in a quantum field theory. We also give full details on the required analytical continuations, 
and we include finite quark chemical potential.^ In particular we investigate quantitatively the effect of Friedel 
oscillations in the screening of static charges at low temperature and ultrarelativistic density. 



* Electronic address: 'hlaizot@Rc t.itl 
^Electronic address: ipp@ect.it 
Electronic address: rebhana@hep.itp.tuwien.ac.at" 
An extension of the large-A^j^ gluon self-energy to finite mass, albeit in the weak coupling regime, can be found in Ref. where 
transport coefficients are calculated in the large- A''/ limit. 
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Figure 1: In the large- A^/ limit, the leading contribution (~ A^^) to the bosonic self-energy H^^ is given by this diagram with 
one fermion loop insertion. 



II. LARGE-iV/ QED AND QCD 

In the large- TV/ limit, an effective coupling g^^j = (for QCD; g^g = e^Nf for QED) is introduced which is 

kept of the order g^^ ~ as TV/ ^ oo and g^ 0. All gluon-gluon interactions of QCD are suppressed by inverse 
powers of Nf - the theory becomes QED like. Also quark-gluon (or electron-photon) interactions are suppressed 
unless they form a new fermion loop. The leading contribution to the gauge boson self energy is therefore simply 
given a fermion loop (Fig-QJ, which is trivially gauge independent. Dyson resummation of this contribution turns the 
gauge boson propagator into a nonperturbative object. 

At large TV/, QCD loses asymptotic freedom. The scale dependence of the coupling in this limit is completely 
determined by the one-loop beta function 

P{g,s) = -g^ - ■ (1) 

Further corrections are suppressed by at least one power of 1/TV/. Integrating this differential equation yields 



(2) 



which impHes a Landau pole at A^, ~ fie^^^ /daitifi) ^ 

The presence of a Landau pole means that large- A^/ QED and QCD only exist as a cutoff theory. However, at finite 
temperature and chemical potential, where one is primarily interested in thermal effects, the presence of a cutoff can 
be neglected as long as T, /i ^ A^. This requires that g^g with renormalization scale of the order of temperature T 
or chemical potential /x, respectively, has to remain smaller than ~ 36, giving enough room for testing the extension 
of hard-thermal-loop expansions from weak to stronger coupling. ^ 



III. LARGE-iV/ GAUGE BOSON PROPAGATOR 

For definiteness, we shall use the Coulomb gauge and decompose the gauge propagator into a longitudinal and a 
transverse contribution 

Goo(Q) = Gl(Q), 

G.,(Q) = |<5.,-^|gt(Q), (3) 
where the corresponding self-energy components are defined through 

q-^ + TIl(Q) 



^ The consequences of a Landau pole in thermal scalar field theory have been studied previously in Refs. fT^ll9l . 
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with Minkowski metric (+, — , — , — ) for the 4-dimensional momenta Q = {qo, q), i.e. = q^ — q^. The corresponding 
spectral functions are determined through the discontinuity of the Feynman propagator along the Minkowski axis 

p{qo,q) = - iG{qo + ie, q) - G{qo - ie, q)) . (5) 

Propagating modes correspond to peaks in the spectral function. In the case of HTL, those peaks are given by infinitely 
narrow delta functions, but beyond the HTL approximation the peaks will acquire a finite width. Such peaks usually^ 
indicate the vicinity of a pole of the propagator to the real axis, but since the propagator does not exhibit a pole in 
the physical sheet, we have to extend our search to the unphysical sheet by analytic continuation. The position of this 
pole can be extracted numerically from the zero of the inverse propagator by specifying either q and extracting qo (q) 
from G~^{qo{q), g) = or by specifying qo and solving G~^{qo, q{qo)) ~ 0. In the former case, we obtain information 
about the frequency uj{q) — Keqo{q) and the decay constant j{q) = linqo{q) for a propagating mode with given 
real-valued wave- vector q, which directly appears as a peak in the corresponding spectral function. In the latter case, 
we obtain the complex wave-vector of a mode induced by a perturbation with a given real-valued frequency uj. 
In the static limit w ^ we can extract the Debye mass from the longitudinal component Gj^^{qo 0, ztod) = 0. 
The transverse component is not screened in the static case, but exhibits dynamical screening for non-zero frequency. 



A. Analytic continuation and scale dependence 

The main technical difficulty in finding the poles of the propagator of Eq. Q, or correspondingly the zeros of the 
inverse propagator, is the correct analytic continuation of the self-energies from the physical sheet into the neighbouring 
unphysical sheet. The propagator G{uj + ie, q) neither contains poles for e > nor for e < 0, but these two analytic 
regions are separated by a discontinuity along e = 0. In order to obtain the correct pole, one has to analytically 
continue the function G{ijj + ie, q) across this separating line from e > to e < 0. 

The self-energy components that appear in the propagator 1^ can be split into a thermal piece and a vacuum piece 

nL(Q) = nL,th(Q)-^nvac(Q), 

nT(g) - nT,th(Q) + nvac(Q), (6) 

where the thermal pieces can be expressed as one-dimensional integrals (see Appendix El, while the vacuum part 
depends on a renormaHzation scale jl, and in the modified minimal subtraction (MS) scheme is given by 

nvac(g) = 3^Q'|iog(^)-^^ (7) 

The analytic continuation of this vacuum piece is straightforward, as one just has to move the branch cut of the 
logarithm appropriately. One possibility is to split up the logarithm log(— Q^) log(g -|- qo) + \-Og[q — qo). As the 
logarithm has a branch cut along the negative real axis, we would cross it for Re g±w < when continuing go = oj + ie 
from e > to e < 0. In the first logarithm we should therefore place the branch cut along the negative imaginary 
axis, and in the second logarithm along the positive imaginary axis. It is useful to define 

logT(z) log(iz)- y, (8) 

logJ'(z) log(-iz)-Hy, (9) 

where log(z) is the standard logarithmic function with a branch cut along the negative real axis. Then log^(2;) has 
the branch cut along the positive imaginary axis, log'''(z) has the branch cut along the negative imaginary axis, while 
they both agree to log(z) = log^(z) = log^(2;) for Re(z) > 0. The analytically continued vacuum contribution can 
then be written by replacing 





log^ [q + qo) + log'f {q - qo) - log {fj^) (10) 



^ See however Ref. I20I1 for an interesting counterexample. 
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in equation 0. 

The analytic continuation of the thermal self-energy functions turns out to be more involved. The expressions for 
Hx.th and IlL.th, which are given in Appendix El involve a one-dimensional integral over a real variable k that has to 
be evaluated numerically. One might assume that it suffices to rewrite the logarithms appearing in the integrands, 
similar to the vacuum contribution Ilvac, and leave the integration over the real variable k untouched. It turns out 
that this prescription does not give the correct analytic continuation for the functions IlT,th{qo,q) and nL,th(9o,9), 
as one can convince oneself by observing a discontinuity in the second derivative of the self-energies at e = 0. Indeed, 
we have to be careful that the integration path for k stays away from logarithmic singularities or branch cuts as we 
approach and penetrate the border of the physical sheet. In order to obtain the proper analytic continuation we 
have to deform the /c-integration path such that crossings with singular points from the logarithms will be avoided. 
Another subtlety involved is that the logarithmic arguments first have to be rewritten in such a way that it is possible 
at all to deform the integration path satisfactorily. Details of this procedure are given in Appendix IbI 

For later use we note that the general form of the analytically continued inverse transverse propagator can be 
written as 

Gt\Q) = -Q' + <?cffrV(|) + nvac(g) (11) 

where the renormalization scale /2f^ appears explicitly in the temperature-independent vacuum piece Ilvac and im- 
plicitly in the coupling g^g — 5off(/^Ms)' scalar function f{Q/T) — IlT{gcS,T,Q)/{g^gT^) is independent of 
the renormalization scale ftj^ as 11^ is proportional to g^g in the large- A^/ limit. 

In the hard thermal loop (HTL) or hard dense loop (HDL) limit H^l, Hl and IIt are given by elementary 
functions (Hsted in Annendix lA 2ll . The HTL limit is derived by assuming soft external lines Q ^ gT . In this 
kinematical regime, the leading contribution is not only gauge independent (also beyond the large- A^/ limit), but also 
formally independent of the renormalization scale, as the vacuum contribution Hvac is suppressed by g^ compared to 
the leading HTL result. This leading piece is of the order of the Debye mass (HTL/HDL quantities will be marked 
by a hat) 

which is the only scale in the HTL/HDL propagator. The HTL plasma frequency, the frequency above which there 
exist propagating longitudinal and transverse modes, is given by a)pj = m|,/3; the effective mass of transverse modes 
in the limit of large momenta (asymptotic thermal mass) by to^ — m|)/2. 

In the \&vge-Nf limit, there is no such simple scaling, and the shape of the dispersion relation could in principle 
depend on the coupling, the temperature, and the renormalization scale. We have a restriction through the exact 
renormalization dependence given by Eq. (QJ which means that we can basically choose a scale for plotting the 
dispersion relations as the result does not change if we vary renormalization scale and coupling according to the 
/3-function. 

To see how this works, let us assume that we know the position of a pole of G{Q), i.e. G^^{Q) — 0. We want to 
show that fLdG~^{Q)/dfi = is also true, provided that gcS and fL follow the /^-function Pig^g) = p.dg^ff/dfi. We 
find that the derivative of ltTT)l 

^ d-^ = '^ + ^(5cff)(r/(-) + ^(ln^--) )+^^(-2) 




/3(5eff) +^ -^=0 (13) 



vanishes indeed, where we used G^^{Q) = to get from the first to the second Hue, and the large- A^/ /3-function 
Eq. in the last Hne. Stated the other way round, we could have derived the /3-function by demanding that 
dG~^{Q)/dp, = vanishes. A similar proof holds of course for the longitudinal component as well. Also the solution 
uj{q) given implicitly by G~^{LLj(q), q) = is independent of ft^jg- Since the Landau pole is also a solution to the zeros 
of the inverse propagator G^^(ALandau) = 0, the position of the Landau pole also obeys the renormaHzation group 
equation. The same is also true for the Debye mass, which is the solution at zero frequency oj and purely imaginary 
q for the longitudinal propagator G^^iuJ — 0,q = irriD) — 0. 

Therefore, if we plot the large- A'/ dispersion relations for a specific choice of the renormalization scale, we actually 
plot the curve for all couplings and renormaHzation scales that are connected to each other by the renormalization 
group equation. The only problem that we might run into is the presence of the Landau pole. 
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Figure 2: Integration path for the spectral function with three poles: One of those originates from the denominator l/{qo—uj), 
while the other two poles are the unavoidable Landau poles. 



B. Sum rules including the Landau pole 

Because of the presence of the Landau pole in the large- A^/ theory, the basic spectral function sum rule 

G(c.,,)=r^^, (14) 

vaHd for complex cu, has to be modified. In the following, we will calculate the missing contribution. Let us first 
look at the zero temperature gauge field propagator. At zero temperature and zero chemical potential, Uriqa, q) and 
nL(<ZOj q) vanish, and one can analytically calculate the position of the Landau pole of the propagator I^J using only 
the vacuum self-energy Solving for G^^{qQ, q) = we obtain 

9 =Aiiisexp(^-^ + -j = Ai. (15) 

We define the spectral function as the discontinuity of the boson propagator along the real 90 axis 

PT{qo,q) = -{G{qn+ie,q)-G{qo~ie,q)) 
I 

= -{G{qo+ie,q)~G*{qo + ie,q)) 
I 

= 2ImG{qo + ie,q). (16) 

The integral over the spectral function 

dqop{qa,q)_ dqo f G{qo + ie,q) G{qo-ie,q) 



27r qo~Lj 7_oo 27r V qo - uj qa - 



(17) 



can be calculated by closing the contour by a harmless grand half circle at infinity and picking up all pole contributions. 
The Landau pole Ijl5|l provides additional poles in the go plane 



go = -q' = ±«Ai(q). (18) 

The first pole lies in the upper half of the complex plane and its residue can be calculated as: 

G{qo + ie, q) _ Gn^ 1 1_ 

^eb,o^.A,(,) _ ^ - ^^^^^^ _ ^ ^^^^^^ . (19J 
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The pole in the lower half qo = —iAL{q) will receive an additional minus sign from the integration orientation, so that 
the total contribution (including the usual pole) is given by 

dqop{qo,q) n( ^ ■ \ ^^^^ ^ ^o^^ 

G{uj + ie,q)- (20) 



, 27r go - w ' gl^ A| - q^ 

where we have substituted back A|(g) = A| — with as defined in lfT5|l . 

With a full one- loop self energy insertion, the position of the Landau pole l^hiq) as the solution to G^^(go — 
AL{q),q) = for a given q can in general only be calculated numerically. In the previous section we simply had 
Al((j') = a/ A|^ — , but the dependence on q is more complicated in the general case. Following the previous 
derivation for a general propagator, we can write the residue at the Landau pole as 



G{qo + ie,q) 1 ( dG ^{qo,q) 



(21) 



go=«Ai,(«)/ 



Using the pole in the lower half plane and the property that the derivative of G{qo,q) is antisymmetric along the 
Euclidean axis dG~^ {qo, q) /dqo\^^^_^j^^^^^= — dG~^ {qo, q)/dqo\^^^^j^^^^^ {G{qo,q) is real and symmetric on the Eu- 
clidean axis, so its derivative is antisymmetric), the full relation is then given by 



dqa p(go,g) _ , . . 2AL(g) / dG ^{qo,q) 



(22) 



qa — iALiq) y 



where the last part is real, despite the appearance of expHcit "?"s. The last term corrects the usual spectral function 
sum rule due to the presence of the Landau pole. A rough estimate of the size of this correction can be obtained from 
the last term of (EOJ: For small Q'^ it is around 5 x 10~''T~^ for ffofflMMS = = ^ around 2 x 10~^r~^ at 
g^g{TrT) = 36, and increases with increasing Q^. For soft momenta Q ~ gcsT one could estimate roughly that this 
term gives a correction to the propagator of the order of (T/A^)^. For hard momenta Q ^ T the correction is of 
order (T/geff Al)^. 



IV. RESULTS 



A. Dispersion relations at zero chemical potential 

Figure El recalls the pressure obtained in the large- Nf limit for finite temperature and zero quark chemical potential 
The full line gives the result at large Nf. The ambiguity introduced by varying the numerical cutoff 
below the Landau pole Q,^„ax = '^^l between a = 1/4 and 1/2 is shown as a tiny (red) band for large coupHngs. 
This ambiguity is suppressed by a factor {T/Al)^. Practically, this means that we can study a range of couplings 
desiP'MS ~ ^ safely, where this number of course depends on the choice of renormalization scale p.y[g [a 
change of this scale changes the coupling according to Eq. J^J]. Perturbation theory ceases to work much sooner and 
shows large scale dependences for couplings already at g'^g > 4. Using an optimized renormalization scale, namely 
fastest apparent convergence of the effective mass parameter ttie (FAC-m) , the perturbative result to order g^g follows 
the exact result up to g'^g « 8. Still, there is a large range of couplings 8 < ffoffC"'^) ^ to explore where strict 
perturbation theory fails, but the ambiguity introduced by the vicinity of the Landau pole is still negligible. 

One of the remarkable and so far not easily understood features of the pressure curve is the minimum appearing 
at a coupling g1g{T:T) « 14. Paradoxically, in this strong coupling range g^g > 14 the pressure seems to approach 
again the interaction-free pressure value. Since the nontrivial content of the large-TV/ pressure can be expressed 
entirely in terms of the gauge boson propagator 0, 0, Q| ; it is natural to ask whether this nonmonotonic behavior 
is associated with qualitative changes in the properties of the gauge boson propagator. We shall find that this is not 
really the case. As will be shown in a subsequent paper the nonmonotonic behavior is instead associated with 
nonmonotonic behavior of the asymptotic quark mass, which is a required ingredient when expressing the entropy in 
terms of quasiparticle spectral data ;23i..24i|. 

In the following we will plot spectral functions and dispersion relations covering weak and strong coupling for 
the gauge boson propagator in the large-A^/ limit, and compare them to the corresponding HTL results. For larger 
coupHngs, the question of renormaHzation scale dependence becomes important. In principle one could plot the 
following comparisons between exact large-A^/ and HTL results at various scales. To avoid cluttering the plots, we 
decided to scale large- A^/ quantities by the large- A^/ Debye mass mu, and HTL quantities by the HTL Debye mass 
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Figure 3: Exact result for the interaction pressure at finite temperature and zero chemical potential at large Nf as a function of 
5off(/^MS" ~ ''^T), compared to strict perturbation theory 12, 13, 14]. The tiny band appearing for large values of the coupling 
for the exact result shows the cutoff dependence from varying the upper numerical integration cutoff between l/\/4 =50% 
and l/\/2 ^70%, of the Landau pole Al. The results of strict perturbation theory are given through order g^g (dotted line), 
QcB (dashed), g^a (dash-dotted), and geff (dash-dot-dotted) where the renormalization scale fly[g is varied between ^tvT (line 
pattern slightly compressed), ttT, and 2ttT (line pattern slightly stretched). The line labelled "FAC-m" indicates the scale 
chosen by the prescription of fastest apparent convergence as indicated in the text for which the curves of g^g and g^g coincide. 



mD/rriD niD/T 

0.5 1 1.5 2 2.5 3 3.5 




Figure 4: Comparison of the large- A^/ Debye mass mo to the corresponding HTL value mo- The renormalization scale /i^jg is 
varied by a factor of 4 around FAC-m. Using FAC-m or ttT as the renormalization scale brings the HTL result within a few 
percent of the full result for the whole range of couplings g^ft plotted. Note that the abscissa does not denote the coupling at 
a fixed scale, but at a different scales as indicated on the lines. 

iriD and present them in the same plot. Each of the quantities is then separately renormalization scale independent, 
while the renormalization scale dependence of the corresponding HTL approximations can be read off separately from 
a plot TOri/m£). (Note that the plasma frequency would be less appropriate to set a scale of the order gT because 
of the imaginary part it obtains at larger couplings g^g - the Debye mass on the contrary always stays purely real.) 
For the spectral functions we note that g'^ffG{qo,q) and therefore glsPiqa^q) is a renormalization scale independent 
quantity. 
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Figure 5: Real and imaginary part of the plasma frequency ujpi in the large- A^/ limit compared to the HTL value Cjpi/rriD = 
l/\/3 ~ 0.577 as a function of 5cff (Mms' — ""2^) at the renormalization scale ttT. The imaginary part of the plasma frequency, the 
damping constant, is caused by the vacuum process of quark-antiquark pair creation, and thus missing in the purely thermal 
HTL approximation. 



One of the most prominent properties of the gauge boson propagator is the emergence of a screening mass in the 
static longitudinal sector, the Debye mass. Numerically, we obtain the Debye mass as the solution to Eq. I^J in the 
static limit, 



9^ -I- nL(go = 0, g) = aXq = imD- 



(23) 



While in the HTL approximation nL('i'o = 0, g) = ?fi|, is a constant given by Eq. Ijl2|l . this is no longer the case for 
the complete expression in the large- A^/ limit. Nevertheless, the solution for mo indeed turns out to be a purely 
real number. The latter can be identified with the screening mass in a Yukawa potential for static sources. Figure 
shows a comparison of the Debye mass in the large- A^/ limit, normaHzed to its corresponding HTL value lfT2|l . Other 
than the following plots, the axis of Fig. 01 is not labelled at a fixed renormalization scale glglirT), but rather at a 
different scale for each curve as denoted next to each line. The four curves are in fact related by the renormalization 
group equation IpJ, because the position of the large- A^/ Debye mass is renormalization scale independent, like any 
singularity of the propagator. The HTL result on the other hand is a truncated result at order g^g. Therefore the 
HTL Debye mass does not run properly with the renormalization scale, leading to a renormalization-scale dependence 
of niD/mD- 

For small an optimal choice of renormalization scale is obtained by requiring that the g^g contribution to the 
Debye mass uid for the full large- iV f result vanishes. This prescription of fastest apparent convergence will be labelled 
as FAC-m and leads to 



-FAC-i 



g5-7B7rr w 0.92 577rT, 



(24) 



which is very close to the choice jly^ = ttT, the lowest Matsubara frequency. Even for the largest couplings ~ 36 
the deviation from the HTL Debye mass is only of the order of 5% for this choice. It will therefore suffice to compare 
all following quantities at a fixed scale /2j^ = ttT, keeping in mind that a change of the scale e.g. by a factor of 2 will 
introduce additional deviations as displayed in Fig. ^ 

Another quantity of central importance is the plasma frequency: It indicates the lowest frequency for which there 
exist propagating modes in the plasma. Transverse and longitudinal modes share the same plasma frequency in an 
isotropic plasma. Frequencies below the plasma frequency are associated with spatial screening. The plasma frequency 
and the frequency of propagating modes are obtained from 



(25) 



While the dispersion relations obtained in the HTL approximation involve only real frequencies i^{q), the full large-iVj 
result is complex. Real and imaginary parts are plotted separately in Fig. [3 Compared to the real part, the magnitude 
of the imaginary part remains small for all values of the coupling. For small couplings, the plasma frequency approaches 
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Figure 6: Real and imaginary part of the dispersion relation ^{q) as a function of real q for g^^ = 9 in the large-A^/ limit 
compared to the HTL result. Abscissa and ordinate for the exact large Nf result are scaled by the large- A^/ Debye mass mo, 
while the HTL curve is on both axis scaled by the HTL mass mjj. At this moderate value of the coupling, which lies at the 
border of applicability of strict perturbation theory, the plasma frequency is slightly lowered and obtains a small imaginary 
part. The light gray line indicates the lightcone. 
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Figure 7: Same as Figure |21 for the larger coupling g1g{iTT) — 36. Real and imaginary part are further lowered compared to 
the case of (j^ff (tT) — 9. As the longitudinal branch approaches the lightcone faster, also its imaginary part approaches at 
smaller values of q/mo- 



the HTL value, which is given by ujpx. — rhBl^/i ~ 0.57735 mu. As the coupHng increases, the renormaUzation-scale 
independent ratio of plasma frequency to Debye mass decreases by 10% at 5off(MMS ~ ^'^) ~ ^^^^ 
g^g — 36. The imaginary part is due to the vacuum process of quark-antiquark pair creation, which is missing in 
the HTL approximation. (In QCD at finite Nf, the damping rate of propagating bosonic modes is dominated by 
Bose-enhanced gluonic contributions proportional to a single power of g rather than to g^g = g'^Np/2 HEEEESl-) 

Figures El and 13 show the dispersion curve for representative values of the coupling glg(jl = ttT) = 9 and 36. These 
plots are scaled by the large- iV/ plasma frequency and thus show a renormalization scale independent result. HTL 
quantities are scaled by the HTL Debye mass. One notes the lower plasma frequency compared to the corresponding 
HTL value, albeit the effect is moderate even at large coupling. The imaginary part of the longitudinal branch vanishes 
quickly as the curve approaches the light-cone. In the limit of vanishing q/iriD, the real as well as the imaginary parts 
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Figure 8: Real and imaginary part of the dispersion relation q{ui) as a function of real lj for g^g — 9 and 36 (labeled in the plot 
as 3 and 6 respectively for Qcs) in the large- A'^/- limit compared to the HTL result. The left panel shows the dispersion relation 
for the transverse component qriu;), while the right panel shows the longitudinal contribution The exact laxge-Nf result 

is scaled by the large- A^/ Debye mass mn while the HTL result is scaled by the HTL Debye mass rfiD. Notably, the real part 
of q{(jj) vanishes in the static case a; in the laxge-Nf limit also for larger couplings g^g, and one can see from the right 
panel that the Debye mass stays a real quantity even beyond HTL. 



of transverse and longitudinal contributions coincide.^ The unique plasma frequency at g = evolves according to 
Fig. 13 

Instead of assuming real wa.ve vectors q and solving numerically for complex uj{q) as we had done in Fig. El we could 
start from a given real frequency to and calculate the corresponding complex wave vector This corresponds to 

calculating the Hnear response to perturbations with a given real frequency. The result is shown in Figure |H1 for two 
different coupHngs 5off(MMs) = ^ (which are labelled as ^off = 3 and 6 respectively in the plot). Above the 

plasma frequency, one has weakly damped spatial oscillations, and below the plasma frequency there is predominantly 
exponential screening, modulated by comparatively long-range oscillations that are absent in the HTL case (though 
fundamentally different from the Friedel oscillations [23; 1^ we shall discuss further below) . Following the curves 
down to vanishing frequency lu ^ leads to the (real) value of the Debye mass in the longitudinal mode (right panel 
of Fig. 13, and in the transverse mode (left panel) to unscreened static modes (but significant dynamical screening 
at small but non- vanishing uj). The real part of (?(w), which corresponds to spatial oscillations, vanishes for both 
transverse and longitudinal modes, in the static Hmit w ^ 0. The imaginary part of the longitudinal modes for 
gcs = 3, 6, and HTL, approaches the value 1 in the static case as the curves are normaHzed to this value. The HTL 
results of Fig.|Hl which are given by dotted and dash-dotted lines, are of course independent of the coupling in units 
of the HTL Debye mass. In a proper comparison between large Nf and HTL one should keep the renormalization 
scale dependence oiniD/mD from Fig. 31 in mind. 

Figures El and El show the spectral functions at two different values of q/mo- The combination T'^g1^p{ijj, q) on the 
ordinate is chosen in order to show a renormalization scale independent result: The quantity g^f^G is renormalization 
scale independent, which can be seen by combining Equations JSJ and ljll|l . In all four plots the broadening of the 
propagating mode peak with increasing coupling is evident. For the smaller wave vector q/mo = 0.2 in Fig. ^1 one 
furthermore nicely observes the decrease of the plasma frequency with increasing coupHng. The HTL pole corresponds 
to a delta function peak, indicated by the vertical line. 



■* This is due to the isotropy of the plasma. See Ref. |3 for counterexamples in anisotropic plasmas within the hard-loop approximation. 
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Figure 9: Spectral function pT and pl as a function oi Lo/niD for q/mo = 1 for various couplings ffcff (/^ms" = '""2^) (^^i^ lines in 
the plot are labelled by the values for (?cff ~ 2, 4, 6). As mentioned in the text, g'^f{p{Lo, q) is a renormalization scale independent 
quantity, and T'^gls scales like rh\, such that the HTL results for various couplings lie on top of each other. The position of the 
infinitely narrow HTL 5-peak is indicated by the vertical line, while the peaks at large Nf have a finite width and are shifted 
to lower energies with increasing coupling. 
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Figure 10: Same as figure lO for for q/mo ~ 0.2. Due to the imaginary part of the plasma frequency, the peak of the 
propagating modes broadens. Additionally their frequency lowers. 



B. Dispersion relations at finite chemical potential 

In the following, we extend the analysis to finite chemical potential /i and zero temperature. Fierure ITTl shows the 
corresponding Debye mass. Here the prescription of fastest apparent convergence gives 

FAC-m = 2Ve/i ~ 3.2974/x at T 0. 

The FAC-m scale for arbitrary finite T and ^ can be found in Ref. In addition to multiples of the FAC-m scale, 
also the simple scale choice = 2fi, varied within a factor of 4, is depicted. 

For the propagating modes, Fierure [121 shows the imaginary part of the dispersion relations at zero temperature for 
the two couplings g^g = 9 and 36 at the scale = 4/i. We did not plot the real part of the function separately, 
since they are almost indistinguishable from Figures El and [3 at zero chemical potential. The HDL result for the 
imaginary part in Fig. El is just zero, and the resulting imaginary part beyond HDL shows a threshold around the 
Fermi momentum g > ^ which can be explained as follows. 

The process responsible for the imaginary part of the self-energy and thus for the propagator pole is the pair creation 
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Figure 11: Comparison of the laxge-Nf Debye mass mo to the corresponding HDL value rho at finite chemical potential and 
zero temperature. The FAC-m scale takes the value fi-fjg = 2^jj,, and is varied by a factor of 4 in the plot. For comparison, 
also the scales /i^jg- = and 4/i are displayed. The full Debye mass is in remarkable agreement with the HDL result for 

the FAC-m scale. 
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Figure 12: Imaginary part of the dispersion relation Lu{q) at T = as a function of real q for (;eff(MMS ~ ^m) — 9 and 36 in the 
large- A^/ limit compared to the HDL result. In contrast to the imaginary part at zero chemical potential in Figures Hand El 
the imaginary part is zero up to a threshold q > 2fj, — Rea;(g) as explained in the text. The constant HDL line in the right 
panel is hidden by the exact large- A^/ lines. The real part of the dispersion relations looks almost identical to the case of zero 
chemical potential (Figures El and and is not plotted separately. 



as shown in FigurelTsl At real momenta q we are looking at time-like momenta for which the Landau damping process 
is kinematically forbidden, even for massive quarks. The only process that could contribute to the imaginary part 
of the self-energy is the creation of a quark with momentum p -|- q and an anti-quark with momentum p. At finite 
chemical potential, this process is kinematically restricted as the Fermi sea is filled up to the Fermi energy Ep = H 
and the newly produced quark must have an energy above, Cp-i-q > The anti-quark can in principle be produced 
within the Fermi sea, Cp > 0, but the dispersion relation of the gluon uj{q) dictates a minimal energy. From energy 
conservation of the process w = tp + Cp+q we obtain for the angle between p and q 
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Figure 13: Diagrammatic processes (a) for Landau damping and (b) for gluon decay. Only the latter one is allowed kinematically 
for time-like gluons with ui'^ > For non-zero chemical potential and zero temperature, this process is further restricted to 
uj + q > 2(1 as explained in the text. At the threshold, a quark with Fermi energy /i and a coUinear anti-quark with energy 
Ep = (a; — g)/2 are produced. 



For zero mass, requiring a valid angle —1 < 008 6* < 1, we obtain the kinematical restriction uj + q>2p>uj — q. For 
a small finite quark mass < = — 9^ , the lower bound for p yields 

P>1 {^l{Q) - q) (27) 

with j{Q) = \/l — A:W? IQ"^ . From the energy threshold of the quark ep+q > ^ we immediately obtain 

uJ-fiQ) + q > 2(1 , (28) 

where 7(Q) = 1 in the massless case. The momentum thresholds in Figure[T2liust follow this formula q > 2fj. — Reuj{q), 
scaled by the Debye mass mo, where real part of uj{q) can be read off from Figures El and H It also explains why the 
momentum threshold for transverse modes is slightly lower than for the corresponding longitudinal mode, because 
^T{q) > ^L{q) for g > 0. 



C. Friedel oscillations 



For finite chemical potential and comparatively small temperatures, electrostatic screening at sufficiently large 
distances is no longer given by a simple Yukawa- type potential involving the Debye mass mu- As is well known in 
nonrelativistic quantum many-body theory , the sharpness of the Fermi surface leads to an asymptotic behavior of 
the potential that is governed by so-called Friedel oscillations, implying a shell-like structure with alternating screening 
and overscreening. At T = O^he amplitude of these oscillations decays by a power law and not like an exponential. 
As has been shown in Ref. j^fll) a similar modification of Debye screening occurs in high-density QED and QCD. In 
the following we shall investigate this effect quantitatively in the large- iV/ theory. 

The free energy of two (Abelian) static charges at a distance r is given by the Fourier transform of the electrostatic 
gauge boson propagator 

d^k if^.f 1 
(2^)3 fc2 + ni(fco -0,A:) 

The deviation from Coulomb law behavior is described by the screening function 



(Ab 

m 

/rj-ih 1 



~ " ^ loo k^+Ii,{k,=Q,kY 

where we have used the fact that 11^(^0 = 0, k) is an even function of k — V^. In the HTL approximation, where 
nL(fco = 0,k) = mjj, one has simply s{r) = e"™"''. 

However, as can be seen from the analytical results in Appendix IB 41 the full gauge boson propagator at zero 
temperature also involves logarithmic branch cuts which touch the real axis at momenta k = ±2/i. Following Refs. |29l 
HDl, we separate these branch cuts infinitesimally according to fc = ±2/Lt ± ie so that the integration path of Eq. ij^Cijl 
is threaded through the four branch points as shown in Fig. El Ignoring (as implicitly done in Ref. j^fll) the Landau 
poles located at real k = ±Al ^ /i by moving them slightly into the lower half-plane, the screening function receives 
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Figure 14: Analytic structure of the electrostatic propagator at zero temperature and integration paths for the screening 
function 18011 . If the Landau poles at real k = ±Al are moved infinitesimally into the lower half-plane, the integration over the 
real k axis can be deformed to encircle the two branch cuts starting at fc = ±2/i + ie and the Debye mass pole at k = imo ■ 

contributions from the Debye mass pole at fc = imo and two integrations along the cuts starting at k = ±2fi + ie. 
Because of the latter, the electrical screening at sufficiently large distances is dominated by oscillatory behavior with 
wavelength n/fj,. In the ultrarelativistic case the asymptotic behavior has been calculated in Ref. (^^I with the result 



which we have transcribed to the MS scheme that we are using here. Nonzero temperature moves the branch cuts 
away from the real axis, leading to an extra exponential suppression factor [s^] e~^'^'^^ so that they cease to govern the 
asymptotic behavior for T > m£)/(27r). We shall therefore concentrate on the T = case, where Friedel oscillations 
should be maximally important. 

With regard to potential implications to high-density QED or QCD it is of interest to know whether deviations 
from standard Debye screening occur already at moderate distances before the screening function has become neghgibly 
small. This requires a numerical evaluation of the integral in Eq. and the results for g'^g = 9 and g'^g = 16 are 
shown in Fig. ^1 The dotted line shows the contribution from the Debye mass pole with values mu = 0.94269/i 
and niD — 1.24705/i, respectively; the long-dashed lines labelled "Friedel" correspond to the asymptotic formula (jHlll . 
The full large- iV/ result which excludes contributions from the Landau poles by moving them infinitesimally into the 
lower half-plane is given by the Hne labelled "exact". As one can see, significant deviations from Debye screening 
appear only at such large distances that the screening function s(r) is of the order 10"'* or below for g^g = 9 and 
g^g = 16. At these couplings, the Landau pole is at scales of the order of 10^/x and 3 x 10^^, respectively. If instead 
of circumventing the Landau pole we cut off the momentum integration at qAl with a varied between 1/2 and 
this affects the screening function already considerably above the scale where Friedel oscillations occur, as shown by 
the dash-dotted lines in Fig.^J^ In real QCD, the one-loop Dyson-resummed gluon propagator also involves vacuum 
contributions from the gluon loop which gives asymptotic freedom and eliminates the Landau pole.® Because of this, 
the results obtained by simply circumventing the Landau pole in the large- A^/ theory are probably a reasonable model 
of screening in high-density QCD at low temperature. Taken as such, these results suggest that Friedel oscillations 
occur only at distances where the screened potential is already very small. At low but nonzero temperature, the 
importance of Friedel oscillations is further reduced by an extra exponential suppression factor e"^'^^'' until they 
disappear completely for T > mD/{2n). 



' A cutofF independent of r actually would give rapid oscillations at large distances. In Fig. 1151 we have instead pinned the cutofF to the 

last zero of the oscillatory integrand below aAj^. 
^ It also introduces a branch cut lying on the imaginary axis, which is however not stable against higher-order contributions. 




(31) 
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Figure 15: Electrostatic screening in the laxge-Nf theory at T = for effective coupling <?efi(MMS~4M) — 9 and 16. The full line 
labelled "exact" shows the numerical result obtained by circumventing the Landau pole in the upper half-plane (see Fig. I14II : 
the lines labelled "cutoff" correspond to cutoffs at the last zero of the integrand in Il3n|l below |fc| — aA_t with a — 1/2 .. . l/\/2. 
The dotted line labelled "Yukawa" corresponds to an exponential with exact large- A^/ screening mass mD{gcs), the long-dashed 
line labelled "Friedel" represents the asymptotic formula 13111 . 



In ultrarelativistic (massless) QED, where g^g ~ 0.092, Friedel oscillations are much more suppressed than in the 
examples above. Ordinary Debye screening dominates out to distances r ~ 250/x~^, where the screening function 
s < 10"^". However, Friedel oscillations in massive QED are less strongly suppressed - they then come with an 
extra factor of mr in the screening function [2^1^. Indeed, Friedel oscillations have been seen experimentally in the 
nonrelativistic electron gas. 



V. CONCLUSIONS 



In this paper we studied the spectral functions and dispersion relations of the gauge boson propagator in the large- 
Nf limit. This allows us to investigate the properties of this propagator not only in the weak coupling regime around 
g'^g{nT) < 4, but also in the strong coupling regime 4 < g1g{-KT) < 36, where strict perturbation theory fails but the 
influence of the Landau pole is still negligible. For larger values of the coupling, the usual spectral function sum rules 
have to be corrected due to the presence of the Landau pole. 

Numerically we calculated the Debye mass, the plasma frequency, dispersion relations, and spectral functions at 
finite temperature and compared those to the analytically known HTL results. The poles of the propagators at large 
Nf are renormalization scale independent. Depending on their vicinity to the real energy axis, they appear as peaks 
of finite width in the spectral functions. These poles do not lie on the physical sheet of the propagator, so a proper 
analytic continuation of the self-energies into the neighbouring unphysical sheet is necessary. In addition to moving 
the branch cuts of the logarithms involved as it is sufficient for the vacuum piece of the self-energy, one also has to 
adjust the integration path of a numerical integration appearing within the thermal piece of the self-energy in order 
to avoid a crossing with logarithmic singularities or branch cuts. 

In the large- TV/ Hmit, the full Debye mass is remarkably close to the HTL value for renormalization scale « /iT 
even up to extremely large couplings. The plasma frequency, and to lesser degree the frequency of the propagating 
modes, acquire an imaginary part which is caused by the vacuum process of quark-antiquark pair creation neglected in 
the HTL approximation. The peaks in the spectral functions get broader with increasing coupling and their frequency 
relative to the HTL value decreases. Overall, there do not seem to be significant qualitative changes in the gauge boson 
propagator as one increases the coupling from small to rather large values. By contrast, the interaction pressure of the 
large- iV/ theory exhibits nonmonotonic behavior. As will be shown in a subsequent paper the latter can be better 
understood in terms of the entropy. While the nontrivial content of the large- TV/ pressure can be formulated entirely 
in terms of the gauge boson propagator, a quasiparticle analysis of the entropy shows that the bosonic contribution is 
unsurprising, whereas the fermionic contributions involve next-to-leading order results for the asymptotic quark mass 
which accounts for the nonmonotonic behavior at large coupling. 

We have also investigated quantitatively the effect of asymptotic Friedel oscillations in electrostatic screening at 
zero temperature and high chemical potential with negligible fermion mass. While Friedel oscillations are a truly tiny 
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effect in ultrarelativistic (massless) QED, their importance increases with effective coupHng strength. However, the 
potential of static charges remains dominated by Yukawa-Uke Debye screening up to distances where the screening 
function has dropped by several orders of magnitude and where the effect of a cutoff to remove the Landau pole would 
be of comparable or greater importance. 

Acknowledgments 

The numerical integration paths have mostly been drawn with the packages AxODRAW |^ and Jaxodraw [3^ . 

Appendix A: BOSONIC ONE-LOOP SELF-ENERGY 

1. Full one-loop result at zero mass 

The vacuum part of the bosonic self energy in the Feynman gauge (covariant gauge with gauge parameter a = 1) 
for Nf flavors and Nc colors, coupling as = g^/iAi:), and the renormalization scale /ijjg in the modifled minimal 
subtraction (MS) scheme is given by 



I I -jC, 



log 



^nA. (A1) 



3 



The thermal part of the bosonic self energy reads 

nL,th(<Zo,9) = g^NfHf + N,Hb). (A2) 



.9V 9^-90 



nT,th(<Zo,g) - ■^i^-:^^^iNfHf + N,Hb) + iNfGf+N,Gb)j, (A3) 
for Nf flavors and TVcColors, where the fermionic functions Gf and Hf can be written as one-dimensional integrals 



1 



Stt^ Jo 

and similarly the bosonic functions Gf, and Hi, 



= t;Z2 I dknf{k)[2k-^ 2iL Li-qokL2], (A5) 



= A r dknk{k)(Ak~l^^LA, (A6) 



7o -'---V 4 q 
Hb = dk ribik) 2k - ^ ^2 L, - qokL2 ) , (A7) 







with 



f 2k + q-qo \ f 2k-q + qo \ 

\2k-q~qoJ \2k + q + qo J 

f 2k + q~qo \ f -q + qo \ , , f 2k - q + qo \ 

^2 = log — - 21og — + log , A9 

\2k-q-qoJ \ q + qo J \2k + q + qo J 

and the bosonic and fermionic distribution functions at temperature T and chemical potential /i 

"fc^^) = e^/T-r ^^^^^ 

"^(^^ ^ I + l + e(^+M)/T + l ) ■ (^11) 

These expressions are valid for all complex qo in a rotation from Euclidean space go = "i^ to Minkowski space go = i^+i^ 
with e > 0. For the analytic continuation into the region with e < the expressions Li and L2 have to be changed as 
described below. 

Apart from their HTL/HDL content, the bosonic contributions l|A6|l and ijATp are gauge dependent (see Refs. 
for one-loop results in gauges other than Feynman). 
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Figure 16: Branch cut and singularity structure of the initial integrand for Gf, Eq. llA4(l . for complex k. The singular points fci 
to ki correspond to solutions of 2fc ± g ± go = from the logarithmic argument in Li, Eq. <A8|l . which are pairwise connected 
by logarithmic branch cuts. For qo = u> + ie, all four singularities will simultaneously cross the real k axis as e changes from 
positive to negative sign. 



2. HTL and HDL 



The corresponding HTL or HDL expressions can be extracted by demanding that the external momenta qo and 
q be small compared to temperature T and/or chemical potential fi (that is to take the leading contribution in a 
1/fc expansion of the integrand, since the main contribution is expected to come from large loop momenta k). The 
relevant physical scale is then given by the HTL Debye mass (HTL quantities are marked by a "-hat) 



- 2 2 

mo = 9 



2Nc + Nf 
6 



In this limit, the equations llA4|l to ljA7ll reduce to 



(A12) 
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27r2' 


Hf = 
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V 6 
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27r2 
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1 + TT log — - — 

2g \ q + qQ 



log 



-q + qp 
q + qo 



(A13) 
(A14) 

(A15) 
(A16) 



In the zero temperature limit Gb and Hb vanish. 



Appendix B: ANALYTIC CONTINUATION OF THE SELF-ENERGY 
1. Integration in the complex k plane 

The fc— integration in the self-energy expressions for Gf^b and Hf^b in equations IIA4|I to ljA7ll runs over real k only. 
One might assume that for the analytic continuation of those functions it is sufHcient to take care of the logarithms 
appearing in Li and L2 and leave the real k integration untouched. It turns out that this simple approach fails. 
For continuing the functions properly in the complex domain, we should think of the A:— integration as an analytic 
integration in the complex k plane. Li and L2 provide singular points and branch cuts in the complex k plane that 
move around, as e varies. As long as the logarithmic singularities or branch cuts never touch the integration line, the 
resulting integral is an analytic, complex function. We will demonstrate the necessary procedure for the branch cut 
structure of Li of G/, but the following considerations can be readily applied to L2 and therefore to Gb, Hf, and Hb. 
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Figure 17: First we break up the branch cuts according to Eq. 



Figure ITfil shows the singularity and branch cut structure for the integrand of G/, Eq. (jA4p . which is determined 
for complex k by the logarithms in Li, Eq. ijASII . For = uj + ie, the logarithmic singularities lie at 



ki = -{-q + uj + ie), 
k2 = ^{q + uj + ie), 

h = ^{-q-uj-ie) = ~k2, 



(Bl) 



so with e > and real q and oj we find that fci and fc2 lie initially above the real k axis, while k^ and fc4 lie below. 
The branch cut structure shown in Figure El corresponds to the initial formula for Li from Eq. ijASII : ki and fc2 are 
connected by a logarithmic branch cut, as are ^3 and fc4. The problem is evident as we try to analytically extend 
the integration from the region e > to e < 0: All four singularities will pass the real axis at the same time, as do 
the corresponding branch cuts that penetrate each other. How can we choose an integration path that avoids contact 
with these singularities and branch cuts? The solution Hes in the fact that we can move around the branch cuts. 

First we break up the branch cuts from the logarithmic expressions in Li and L2 as shown in Figure fT71 Using 
the logarithmic functions log^(z) and log^(2:) from Eqs. © and ^ with branch cuts along the positive or negative 
imaginary axis, we can rewrite equation ijASII as 



- break— up 



logT (2fc + g - qo) - logT {2k-q-qo) 
- log^ {2k- q + qo) + logH2fc + q + qo), 



(B2) 



which agrees with Li on the real k axis, and similarly for L2 

^break-up ^ logT (2fc + g - g^) - log^ (2fc - g - go) 
-2 (logH-9 + qo) - log^g + qo)) 
+ log^ {2k- q + qo) - log^ {2k + q + qo). 



(B3) 



The branch cuts in the two logarithmic expressions of L2 independent of k (second line) are chosen such that the 
functions continue analytically for real q and uj when e changes its sign from positive to negative. Using Li'^'^^ ""'^ 
gives the branch cut structure as shown in Figure El It is now clear how the path should be deformed from e > in 
Fig-E|to e < in Fig. El To further justify this procedure, let us note that one could have used the new integration 
path from Fig. El already in Fig. El where it would not have changed the result of the integration, as the result of 
the complex integration is of course path independent. Keeping the same S-shaped integration path in both figures, 
it is clear that moving a logarithmic singularity across the real k axis is not different from moving it around anywhere 
else in the complex k plane, provided no singularities are crossed. 
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Figure 18: In this way, we can deform the integration path when the singularities and branch cuts move across the real k axis. 
Note that we could have used this integration path already in figureEl In that case it is obvious that moving the singularities 
across the real k axis would not be different from moving the singularities around anywhere else in the complex k plane. 
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Figure 19: If Refci > and Refcs < then the path would look like this. In Figures IT?)! to 1181 we had assumed Refcs > and 
Re ki < 0. 



Practically, to change from positive to negative e, we could either literally choose a numerical integration contour 
similar to Figure ITM or we could basically stick to the original integration along the real k axis, taking correcting 
branch cut contributions from the real axis to, say, fc2 and ki into account. At first sight, the latter strategy seems 
insofar appealing, as the branch-cut contributions can be calculated analytically without the necessity of numerical 
integration. This strategy is briefly discussed below in Annendix IB 31 But since the integration parts along the real 
k axis have to be integrated numerically anyway, it turns out to be easier to choose an appropriate integration path 
for numerical integration in the complex k plane which avoids all branch cuts in a safe distance. 

2. Choosing a complex integration path 

We shall show how to choose an integration path in the complex k plane that circumvents the branch cuts. In 
choosing the path, we only need to care about logarithmic singularities h with positive real part as only those will 
cross the initial integration range k G [0, oo). Assuming that Reg > and real uj > in qo = uj + ie, we see that in 
ijBlll the singularity Re fc2 > always lies in the right half of the complex plane while Re fc4 < lies on the left half. 
Also ki and fcs lie on opposite sides of the imaginary fc-axis, so there are two distinct cases for the path topology: 
Either Refci > as in flgure^l or Refcs > as in flgure^l In the latter case, Refca is always smaller than Rek2, 
but in the former case Re ki could be either smaller or larger than Re ^2 . 

At finite chemical potential the fermionic distribution function n/(fc) from Eq. llAlip provides additional poles at 
k — ±fi + {2n + l)i7rT (n £ Z). Figure EDI shows the particular case < Refci < n < Refc2, but the Matsubara 
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Figure 20: At finite chemical potential /i, the fermionic distribution function ijAllll provides additional poles at the Matsubara 
frequencies k = ±/i + (2n + l)iTTT (n € Z) which have to be avoided. 

frequencies can also lie to the left of Re fci or to the right of Re ^2 . As the initial path for the /c-integration (starting 
from e > 0) lies along the real k axis between the poles k = ^± inT, also after the analytic continuation from e > 
to e < the integration path has to stay between these two lowest Matsubara frequency poles. Similarly, one has to 
take care of the Matsubara frequency poles in the case Rek^ > as in figure El 

There is a lot of freedom in choosing an integration path that fulfills all of these constraints and correctly avoids 
the problematic regions in the complex plane. Our implementation of the complex integration path follows essentially 
those suggested in Figures lTsllTol andEOl where we vary the minimal distances to the branch cuts and singularities for 
testing purposes. One can readily check that the final result is independent of the exact implementation of a correct 
integration path. 

The same integration path can be used for calculating the derivative with respect to qo, by forming the derivative of 
the integrand first. The derivative of the self-energy functions can be used in specifying the Jacobian for a numerical 
root search. The search for the Debye mass requires complex q. Assuming Reg > and Img > this implies 
straightforward modifications to the integration paths: As the real part of q{qo) vanishes in the limit qo ^ 0, basically 
ki and ^2 from Eqs. ljBl|l come to lie vertically aligned on top of each other in Fig. ^1 instead of horizontally aligned 
side by side. 

3. Alternative ways of calculating the analytic continuation 

One possibility (which we did not follow in our numerical evaluation) is to analytically calculate the contribution 
caused by a branch cut that crosses our original integration line. Integrating a function f{z)\og{z) (where f{z) is 
analytic) around a part of the branch cut is equivalent to integrating /(z)27ri along the path of the branch cut, e.g. 
with a < 6 < we can write 



pb—ie 


/'a-\-ie 


•/ a—ie 


Jb+ie 



f{z)\og{z)dz = f{z)2Tridz, 

b 



(B4) 



as the real parts of log(z) = log |z| + ? arg z along the two contour parts cancel each other, while the net contribution 
comes from integrating zLin along the contour. 

For we basically have to integrate nf^b{k) x polynomial(fc) x 27ri. This can be done analytically. Here are some 
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examples for rib 

J nb{k)dk = Tlog (e'^/'^ - l) - fc, (B5) 

/J 2 

knb{k)dk = T^U2 (e*^/^) + fcTlog (l - e*^/^) - (B6) 



/ 



k^nb{k)dk = -2T^U:i (e*-'/^) + 2kT^U2 (e^'^) 



+fc2riog(l-e'=/^)-^, (B7) 



where Li„(z) is the polylogarithm function. Similar expressions can be derived for n/. These functions have branch 
cuts for Re(fc) > which would have to be corrected, e.g. 



j nb{k)dk = Tlog {e^''^ - - fc - 27ri 



Im(fc) 1 
2tiT ~ 2 



(B8) 



where the ceiling function \x\ means the smallest integer > x. Similar adjustments would be necessary for the 
other relations above. Since these adjustments are error-prone, and we have to integrate a large part of the result 
numerically anyway, it turns out to be more straightforward to choose a complex integration path and integrate 
everything numerically, as described in the previous section. 

Another possibility to obtain the analytically continued propagator is to exploit the relation between propagator 
and spectral function, Eq. ()14|l . Apart from the fact that this would introduce another numerical integration over go, 
this relation is only reliable for small couplings where the influence of the Landau pole can be neglected. For larger 
coupHngs, the correction term of Eq. lf22|l would have to be taken into account. 

Yet another way to obtain the analytically continued retarded propagator is given by adding a correction term to 
the Feynman propagator (a detailed description of this procedure can be found in Appendix B of Ref. (36] ) 

G^(w - ie, q) = G{uj - ie) + 2i/(w - ie, q) . (B9) 

Along the discontinuity, i.e. for real the correction term coincides with the spectral function from Eq. Q, /(w, q) = 
p{uj,q), but different from the spectral function, the function I{uj — ie,q) has to be an analytic function. By the 
uniqueness theorem of the analytic continuation, one obtains the unique retarded Greens function G^'{lj — ie, q). The 
tricky point of this approach is to flnd the correct analytic continuation of the spectral function p{qo,q). We can 
not readily apply this method to the large-TV/ limit, as the spectral function is, just like the propagator itself, only 
obtained numerically. The proper analytic continuation of the function /(go, would be just as tough as the analytic 
continuation of G{qo, q) itself. 



4. Analytic result at zero temperature and finite chemical potential 
The fermionic distribution function IIAllll can be expanded for small T as 

nfik,T,fM) = ^iei-ik-p)) + e{-ik + p))) 

+T'^ {S'{k - m) + S'{k + m)) + 0(t4) (BIO) 
where S'{k — /i) is to be understood in the sense of distributions 

5'{k~p)f{k)dk = -f\^ji). (Bll) 



In this limit, the integrations in (IA4|I to IIA7|l can be performed anal ytic ally. The analytic continuation is straightfor- 
ward with the tools presented here. The result at zero temperature |22| can be written using 

R{qo,q) = i2fl + qa + q)log{2^i + qo + q)- {qo + q)logiqo + q), (B12) 
S{qo,q) = (2/i + (7o -2g)(2/x + go + 9)'log(2/i + go + 9) (B13) 
-(go - 2q){qo + qf log(go + ?)• 
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We use the following abbreviation to keep a lot of terms in dense notation 

R±{qo, q) = R{qo, q) - R{qo, -q) + R{-qQ, q) - R{-qo, -q) 



(B14) 



and similarly for S±. (just replace -R by S). The functions C'-p and Hy (the upper index "(0)" denotes the T = 
limit) can then be written as 



AO) 



Gf\qo,q) = ^U + '\^Rtiq.,q) 



Hf\qo, q) ^ (sSAi' + - (-2Vgo (log(go + q) - log(go - <z)) + 5±(go, q)) 



(B15) 
(B16) 



In order to correctly analytically continue these formulae into the region of e < 0, one needs to replace all logarithms 
that contain an initially positive (negative) go to log''' (log^) from equations and Defining R^''^^\qo, q) := 
[R{qo, q) ■with log — > log'''^^^] and similarly for S, we can write 



i?±|((Zo, q) = RHqo,q) - RKqo, -q) + RH-qa, q) - RH-qo, -q) 



and the very same for S. The analytically continued functions g'"P and h'"^^ can then be written as 



Gf{qo,q) 



(B17) 

(B18) 
(B19) 

The next term in the small T expansion of equation IIB10|l is easily calculated using |IbTT)| and can be written as 
(the upper index "(2)" denotes the second derivative w.r.t. T at T = 0) 



1 
2^ 



Hf\qo, q) = + - (-2Vgo (W(go + g) - log^go - q)) + S^\{qo, q)) 



Gf\qo,q) = ^ 



Hf\qo,q) - ^[{2fi + qo){\og{2n + q + qo)-\ogi2fi-q + qo)) 

-2qo (log(go + q) - log(go - q)) 

+ {2fi - qo) (log(2/z + g - go) - ^og{2^ - g - go))] . 
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1 



(B20) 



(B21) 



(2) 

As is regular in the whole go complex plane, we only have to take care of the proper analytic continuation of 



H 



(2) 



Hf\qo,q) = ^ [(2/i + go) (W(2/^ + g + go)-log'''(2/x-g + go)) 
-2go {\og^{qo + g) - log^(go - g)) 
+ {2^1 - go) (logT(2^ + g - go) - W(2m ^ 9 - go)) • 
The analytically continued expansion at small T and finite chemical potential is then given as 

G/(go,g) = (go, g)+r2Gf(go,g) + 0(t4), 



(B22) 



Hs{q^,q) = i/f (go,g)+T^i?f (go,g)+0(n. 
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